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We present a new approach to measure the shapes of galaxies, a fundamental task 
in observational astronomy. This approach is based on the decomposition of a 
galaxy image into a series of orthogonal basis functions, or 'shapelets'. Our choice 
of basis functions, namely the Gauss-Hermite series, has a number of remarkable 
properties under distortions, convolutions and noise, which makes them partic- 
ularly well suited for astrophysical applications. In particular, we describe how 
they can be used to measure the shear induced by weak gravitational lensing, with 
the precision required for upcoming surveys. We also show how shapelets can be 
used to reconstruct images from interferometric observations. Other application of 
shapelets, such as image compression, PSF deconvolution, de-projection and the 
study of galaxy morphology, are also briefly discussed. 



1 Introduction 

The measurement of galaxy shapes is a fundamental task in observational 
astronomy. In these proceedings, we present a new method for shape mea- 
surements described in detail in Refregier (2001). It is based on the linear 
decomposition of each galaxy into a series of localised basis functions with dif- 
ferent shapes, which we call 'Shapelets'. As a basis set, we choose the Gauss- 
Hermite series, whose remarkable properties make it ideal for astronomical 
applications. In particular, we summarise the results of Refregier & Bacon 
(2001) who showed how shapelets can be used to measure weak lensing, with 
the precision required for upcoming surveys. We also describe how shapelets 
can be used to reconstruct images from interferometric data, as described by 
Chang & Refregier (2001). 
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We illustrate the method using images from the Hubble Space Telescope 
and from the FIRST radio survey. Finally, we describe several further appli- 
cations of shapelets and discuss their relevance in the context of astronomical 
applications requiring high precision. 



2 Shapelet Method and Properties 

We begin by summarising the formalism of Refregier (2001) for a description 
of galaxies in the Gauss-Hermite basis set. A galaxy with intensity / (x) can 
be decomposed into our basis functions -B n (x; /?) as 

/(x)=]T/ nJ B n (x;/3), (1) 

n 

where x = (x\,X2) and n = (711,712). The 2-dimensional cartesian basis func- 
tions can be written as -B n (x;/3) = B ni {x\; f3)B n2 (x2', 0), in terms of the 1- 
dimensional basis functions 

B n (x;0) = [2Vn!/3p H n (jj e~£, (2) 

where H n (x) is a Hermite polynomial of order n. The parameter (3 is a charac- 
teristic scale, which is typically chosen to be close to the radius of the object. 
These basis functions are the eigenstates of the Quantum Harmonic Oscillator 
(QHO), allowing us to use the formalism developed for this problem. Sim- 
ilar decomposition into basis functions has been independently suggested by 
Bernstein & Jarvis (2001). The first few basis functions are shown in Figure [jj 
Because these basis functions, or 'shapelets', form a complete orthonormal 
set, the coefficients /„ can be found using 

/oo 
d 2 x /(x)5 n (x;/3). (3) 
-00 

Figure [2] shows an example of the decomposition of a galaxy in the Hubble Deep 
Field (Williams et al 1996). The details of the image are fully reconstructed 
if coefficients up to n — 20 are included; this decomposition thus provides an 
excellent and efficient description of galaxy images in practice. 

The chosen basis functions have a number of remarkable properties. Firstly, 
they are invariant under Fourier transforms up to a scaling factor, so that 

B n (k;(3) =i n B n (k;0- 1 ), (4) 

where tildes denote Fourier transforms. As a result, convolutions (which cor- 
respond to products in Fourier space) can be simply expressed using shapelets. 
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Figure 1: First few 2-dimensional Cartesian basis functions B nit „ 2 . The dark and light 
regions correspond to positive and negative values, respectively. 

Let us consider the convolution h(x) = (/ *g)(x) of 2 functions. Each function 
can be decomposed into shapelet coefficients with scales a, (3 and 7, which are 
then related by 

00 

hn = ^ ' CnmlfmOlt (5) 
1X1,1=0 

where the convolution tensor C nm i{a, (3, 7) can be estimated analytically using 
a simple recursion relation (Refregier & Bacon 2001). 

Shapelets also have simple properties under coordinate transformations 
such as translation, dilation, rotation and shear. For instance, let us consider 
the distortion of an image /(x) under the action of a weak shear 7,, as occurs 
in weak gravitational lensing. To first order in the shear, the distorted image 
can be written as 

/'~(l+7^)/, (6) 

where Si is the shear operator. It is easy to show that this operator takes a 
simple form in shapelet space, namely 

S% = i (a\ 2 - a5, 2 - a\ + a%j , S2 = a\a\ - aid 2 , (7) 

where aj and on are the raising and lowering operators in the QHO formalism, 
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Figure 2: Decomposition of a galaxy image found in the HDF. The original 60 X 60 pixel 
HST image (upper left-hand panel) can be compared with the reconstructed images with 
different maximum order n = n\ +U2- The shapelet scale is chosen to be (3 = 4 pixels. The 
lower right-hand panel (n < 20) is virtually indistinguishable from the initial image. 



for each dimension i — 1,2. Similar operators can be constructed for the other 
coordinate transformations. 



3 Measurement of Weak Lensing 

The weak distortions induced by lensing on the images of background galaxies 
provide a direct measure of the distribution of mass in the Universe. This 
weak lensing method is now routinely used to study galaxy clusters, and has 
recently been detected in the field (Wittman et al 2000; Bacon et al 2000; 
Kaiser et al 2000; Maoli et al 2001; van Waerbeke et al 2000, 2001; Rhodes et 
al. 2001). Because the lensing effect is only of a few percent on large scales, 
a precise method for measuring the shear is required. The original methods 
of Bonnet & Mellier (1995) and Kaiser, Squires & Broadhurst (KSB; 1995) 
are not sufficiently accurate and stable for the upcoming weak lensing surveys 
(see Bacon et al 2001, Erben et al 2001). Thus several new methods have 
been proposed (Kuijken 1999; Rhodes, Refregier & Groth 2000; Kaiser 2000; 
Bernstein & Jarvis 2001). 

The remarkable properties of our basis functions make shapelets partic- 
ularly well suited for providing the basis of a new method for measuring the 



4 



shear (Refregier & Bacon 2001). First, the Point-Spread Function (PSF) of 
the instrument can be modeled by decomposing stellar images into shapelets 
and by interpolating the resulting coefficients across an image. Galaxy images 
can then themselves be decomposed into shapelet coefficients, and the analyt- 
ical form of the convolution matrix (see Eq. can be used to deconvolve the 
Point-Spread Function. From the properties of shapelets under shears (eq. JtJ), 
one can then construct a linear estimator for the shear from the (deconvolved) 
coefficients of the galaxies of the form 

~ fn ~ (/n) (8) 

fin — ^ / f \ \°) 
'-'inm \/m/ 

where f' n are the lensed coefficients, the brackets denote an average over the 
unlcnsed galaxy population, and n is even (odd) for the 71(72) component of 
the shear. The tensor Si nm is the matrix representation of the shear operators 
of Equation ([j]). ft is easy to show that these estimators are unbiased, i.e. that 
(7 m) = Tij if the unlensed galaxies are randomly oriented. These estimators 
can then be combined to construct a minimum variance global estimator for 
the shear over a sample of galaxies. 

To test the method, we used the image simulations described in Bacon 
et al. (2001). These were designed to reproduce the observational conditions 
of ground based telescopes such as appropriate throughput, PSF and noise, 
along with the statistics of the galaxy population observed in HST images. An 
artificial shear was applied to the simulated galaxies, and the shapelet method 
was used on the resulting realistic images to recover the applied shear. The 
recovered shear values from the simulations are plotted as a function of the 
output shear in Figure ||. Clearly, the method is unbiased, and is found to be 
robust with PSF shape and size. It thus provides an improvement over the 
KSB method which was shown to have small but significant instabilities and 
biases (Bacon et al 2001; Erben et al. 2001). 

4 Image reconstruction with Interferometers 

Another application of shapelets is the reconstruction of images from interfero- 
metric data, as described in detail in Chang & Refregier (2001). Interferomet- 
ric data are collected in the uv space, the Fourier-transform of the sky surface 
brightness. The observed quantity is the visibility measured for each antenna 
pair at time t and at frequency v and is given by 

V i} (v,t) = [ d 2 l t] e -2,i[ul+vm +w (^uW-l)] (9) 

J V 1 - W 
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Figure 3: Input shear vs recovered shear for a set of 11 simulations. Note the linear rela- 
tionship between input and recovered shear; the slope is found to be 0.97 ± 0.04 for 71 and 
1.00 ±0.04 for 72. 

where /(l, v, t) is the surface brightness of the sky at position 1 = (l,m) with 
respect to the phase center, and A(l, v) is the (frequency-dependent) primary 
beam. 

In analogy with the shapelet decomposition in real space, we wish to de- 
compose sources directly in the uv space. However, the discrete and finite 
sampling of the uv plane prevents a direct linear decomposition of the vis- 
ibilities Vij into shapelet coefficients. Instead, we simultaneously fit for the 
shapelet coefficients of a collection of sources on the uv plane, by using a x 2 
fit: 

X 2 = (d-M f) T CT 1 (d-M f), (10) 

where d = {Vij} is the visibility data vector , M = {V" S } is the theory matrix 
composed of visibilities corresponding to each basis function, and f = {/ ns } 
is the shapelet coefficient vector. The indices s and n stand for the source 
number and shapelet state, respectively. The covariance error matrix of the 
visibilities 

C = cov[d,d]=((d-(d)f(d-<d») (11) 

is provided by the interferometric hardware. The resulting x 2_n t is linear in its 
parameters and can therefore be performed by simple matrix operations. The 
complex effects of bandwidth smearing, time averaging and non-coplanarity of 
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Figure 4: Reconstruction of three sources from one of the FIRST grid pointings. The 
'dirty' images (which are direct inverse Fourier transformation of the data), CLEAN images, 
and shapelet reconstructions are shown from top to bottom, respectively. The images are 
32" across and the resolution is about 5". 4 (FWHM). The dirty and CLEAN images are 
displayed with a 1".8 pixel size, while the shapelet reconstruction images have 1" pixels. For 
the shapelet reconstruction, we used Wiener filtering and smoothing by a Gaussian restoring 
beam with a standard deviation of 2.3". 



the array can be fully corrected for in the process. 

As an example we consider the observing conditions of the FIRST radio 
survey (Becker et al. 1995; White et al. 1997). Using one of the FIRST point- 
ings, as shown in Fig. [|, we find that our method compares well with CLEAN, 
the commonly used method for intcrfcrometric imaging. Our method has the 
advantage of being linear in the shapelet parameters, of fitting all sources si- 
multaneously, and of providing the full covariance matrix of the coefficients, 
which allows us to quantify the errors and cross-talk in image shapes. It is 
well-suited for quantitative shape measurements which require high-precision. 
In particular, combining with the results from the previous section (Refregier 
& Bacon 2001), our results provide an accurate method for measuring weak 
lensing with interferometers. 
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5 Conclusions 



Shapelets provide a new method to measure the shape of galaxies. The chosen 
basis functions have convenient properties under convolutions and distortions. 
This makes shapelets ideally suited for the measurement of weak lensing shear, 
after correcting for the smearing effect of the point spread function. They can 
also be used to reconstruct images from interferometric data. The method is 
very general and is ideally suited for astrophysical applications requiring high 
precision. In particular, it can be used for image compression, modeling and 
deconvolution of the point spread function, and de-projection from 2 dimen- 
sional to 3 dimensional data. Another immediate application to explore is that 
of the measurement and classification galaxy morphology. Shapelets could thus 
help in unifying the different observational and theoretical studies of galaxy 
shapes presented in this volume. 
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